SAXSFit: A program for fitting small-angle x-ray and neutron scattering data 
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SAXSFit is a computer analysis program that has been developed to assist in the fitting of 
small-angle x-ray and neutron scattering spectra primarily from nanoparticles (nanopores). The 
fitting procedure yields the pore or particle size distribution and eta parameter for one or two 
size distributions (which can be log-normal, Schulz, or Gaussian). A power-law and/or constant 
background can also be included. The program is written in Java so as to be stand-alone and 
platform-independent, and is designed to be easy for novices to use, with a user-friendly graphical 
interface. 
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I. INTRODUCTION 

Small-angle x-ray scattering (SAXS) and small-angle 
neutron scattering (SANS) are well-established and 
widely used techniques for studying inhomogeneities on 
length scales from near-atomic scale (1 nm) up to mi- 
crons (1000 nm). Recently, there has been an increasing 
emphasis and importance of nanoscale materials, due to 
the distinct physical and chemical properties inherent in 
these materials [l], This, together with the signifi- 
cant advances in X-ray and neutron sources, has resulted 
in the dramatically increased use of SAXS and SANS 
for characterizing nanoscale materials and self-assembled 
systems. For example, these techniques are used to inves- 
tigate polymer blends, microemulsions, geological mate- 
rials, bones, cements, ceramics and nanoparticles. These 
measurements are often made over a range of length 
scales and in real time during materials processing or 
other reactions such as synthesis. However, there has 
been less progress in SAXS and SANS data analysis, al- 
though some analysis software is available. For example, 
programs based on IGOR Pro primarily for the reduction 
and analysis of SANS and ultra-small-angle neutron scat- 
tering (USANS) are available from NIST ^. PRINSAS 
has been developed for the analysis of SANS, USANS 
and SAXS data for geological samples and other porous 
media j3|. PRIMUS and ATSAS 2.1 are used primarily 
for the analysis of biological macromolecules in solutions, 
but can be used for other systems such as nanoparticles 
and polymers FISH is another SANS and USANS 

fitting program developed at ISIS 0]. The Indra and 
Irena USAXS data reduction and analysis package devel- 
oped at APS Q is also based on IGOR. Both of these 
latter programs offer several advanced features, includ- 
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ing multiple form factor choices and background reduc- 
tion routines. While these programs provide a powerful 
analysis capability, they can be complicated to use and 
some are based on commercial software. This has moti- 
vated the development of a simple, easy to use analysis 
package. 

In this paper, we describe SAXSFit - a program devel- 
oped to fit SAXS and SANS data for systems of parti- 
cles or pores with a distribution of particle (pore) sizes. 
SAXSFit is easy to use and applicable to a wide variety 
of materials systems. The program is most appropriate 
to low concentrations of particles or pores, due to the ap- 
proximations used, but it does account for interparticle 
scattering within the local monodisperse approximation 
. The program is based on Java and is readily portable 
with a user-friendly graphical interface. The emphasis of 
SAXSFit is to provide an easy-to-use analysis package 
primarily for novices, but also of use to experts. 



II. SOFTWARE DESCRIPTION AND USE 

SAXSFit is written in Java (SDK 1.4.2) and provides 
a graphical user interface (Figure [T]) to select and ad- 
just parameters to be used in the fit, change the plotting 
display and range of data to be used, calculating 'initial 
guess' patterns and running the fit. It uses the algo- 
rithms of a Matlab-based program [l3|. The advantage 
of using Java is to provide a stand-alone program which is 
platform-independent, along with having a user-friendly 
graphical interface. SAXSFit is also available as a Win- 
dows executable. 

SAXSFit can read ASCII data files (comma, space, or 
tab delimited), with or without a non- numerical header, 
which consist of two column (g, I{q)) or three-column [q, 
I{q), error bars) data. Any subsequent columns in the 
data file are ignored. Once a data curve is successfully 
imported from an input file, it is plotted in a separate 
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FIG. 1: Screen dump of the SAXSFit graphical user interface and plot window. 



window and the fitting buttons are enabled. The error 
bars are also plotted if the input file contains them. The 
plot can be manipulated by changing the q-range and 
selecting whether it is log-log or linear. Initial guess and 
fitted curves are displayed when they are calculated. The 
g-min and g-max values of the data to be fitted are also 
shown as vertical lines and can be altered by changing 
the appropriate text boxes. Several fitting parameters 
are available, with the option to fit or fix their values. 

Three distributions of particle/pore sizes are available: 
log-normal, Schulz and Gaussian. These use the same 
two fitting parameters, 'particle/pore size' (rg) and 'dis- 
persion' (cr), and are detailed in Section UlI Al The units 
for the pore size are the inverse of the units of the data 
(i.e. A for data in A~^, or nm for data in nm^^). A 
second size distribution can also be included in the fit. It 
has been shown that the choice of distribution function 
does not dramatically affect the final result [ll|, [l3| . A 
constant background and/or power law Aq"^ can also 
be included. 

Advanced options include the ability to change the 
maximum number of iterations and the weighting of the 
data. Again there are several options: a constant weight- 
ing (w = 1), statistical weighting (w = 1/ I{q)), or uncer- 
tainty weighting {w = l/A/(g)^ - only applicable where 
data error bars (A/(g)) have been imported). 

Two output files are produced, consisting of the fit (a 
two-column ASCII file) and a log file (plain text) show- 
ing the values of the parameters at each iteration of the 
fitting process, and the final result including parameter 
uncertainties, reduced and goodness of fit {E? value). 
These are described in more detail in Section [TnEl The 



final parameters are also displayed on the control panel. 

Users must be aware of the assumptions made in the 
modeling of the data, which uses a hard sphere model 
with a local monodisperse approximation. Strongly inter- 
acting systems, for example systems with a high degree of 
periodicity, are outside the scope of this approximation. 
The user is responsible for understanding the applicabil- 
ity of this approximation to their system, and ensuring 
that the fitted results are physically meaningful. 

III. MATHEMATICAL DETAILS 

The small angle scattering intensity I{q) is related to 
the scattering cross section ^{q) by 

I{q)^^,At{m^{q) (1) 

where 0o is the incident flux (number of photons, or neu- 
trons, per area per second), A is the illuminated area on 
the sample, t is the sample thickness and A57 is the solid 
angle subtended by a pixel in the detector p^ . 

For SAXS, the scattering cross-section is calculated 
from the structure factor and particle/pore size distri- 
bution from 

J poo 

— iq) = rliApfN n{r) [f{qr)f S{qr)dr (2) 

where r,, is the electron radius, Ap is the electron den- 
sity contrast, N is the number density, n(r) is the num- 
ber fraction particle/pore size distribution (normalized 
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so that the integral over r is unity), f{qr) is the spheri- 
cal form factor, and S{qr) is the structure factor. These 
terms are defined in the following sections. The final 
equation for the intensity used by the program is 



2. Schulz distribution 



I{q) =c n{r) [f{qr)] S{qr)dr 



(3) 



where the scale factor c is a fitted parameter, equivalent 
to 



(j>oAt{m)rl{^pfN 



(4) 



For SANS an expression similar to Eq. [4] holds. 

The data are modeled using a hard-sphere model 
with local monodisperse approximation [9^. This model 
assumes that the particles are spherical and locally 
monodisperse in size. In other words, the particle po- 
sitions are correlated with their size. This is a good ap- 
proximation for systems with large polydispersity and 
the approximation provides meaningful results providing 
the particles are for the most part not inter-connected 
(e.g., the particle concentration is not too high) and are 
not spatially periodic. For porous systems (with not too 
high pore concentration), the 'particle' radius is equiva- 
lent to the pore size. 



A. Particle/pore size distribution, n(r) 

The user has the choice of three pore/particle size 
distributions, which use the same fitting parameters rg 
('pore size' radius) and a ('dispersion'). If two size dis- 
tributions are selected, the distribution function is ex- 
panded to have the form: 

n{r) = (1 - A)ni(r, ro,i,cri) -f An2(r, ro,2, 0-2) 

where A is the number fraction of the second distribution 
and ro.j and aj are the rg and a parameters for the f '^ 
distribution. For example, a 50:50 mixture by number 
fraction would have A = 0.5. To model a situation in- 
volving a mixture with differing contrasts, A would be 
weighted by the different contrast values. The user has 
a choice of three distributions, as follows. 



Log normal distribution 



n{r) = Z^X 



,_, expi-ZX) 
rT{Z) 



where Z — X — and r(Z) is the Gamma func- 
tion, defined Va; e R by T{x + 1) ^ xT{x), and r(l) = 1 
(14)] . The Schulz distribution is frequently used in SANS 
analysis. It is physically reasonable in that it is skewed 
towards large sizes and has a shape close to a log-normal 
distribution. As Z ^ 00 it approaches a Gaussian distri- 
bution iJj]. 



3. Gaussian distribution 



n(r) 



cr\/27r 



exp 



r - ro 



The Gaussian distribution is symmetric about the 
mean, rg, and has variance (rgtr)^. In practise it is only 
useful for systems with low polydispersity (small a). 



B. Spherical form factor, f{qr) 

The spherical form factor has the following form P.[l6l|: 



f{qr) = rnr^ 



sin(qr) — co8{qr)qr 
(gr)3 



C. Structure factor, S{qr) 

The structure factor follows the local monodisperse ap- 
proximation (LMA) for hard spheres 0, [l3| , given by 



SiqRHs) = [1 + 2ArjGiqRHs)/{qRHsT' 

where rj is the dimensionless parameter eta (sometimes 
referred to as the hard sphere volume fraction, hav- 
ing a value between and 1), Rhs is the hard sphere 
pore/particle radius, defined as Rhs — Cr, where C re- 
lates the hard-sphere radius to the physical particle ra- 
dius 0, and G{qRHs) has the form: 



n(r) = exp -- 



In 



1 



rtrv 27r 



This has a maximum a.t r — tq exp {—cr^ 
rg exp f-) J and variance [exp (2cr^) — exp (ct 



a mean of 

2M 



G{A) = a{sinA- AcosA)/A'^ + 

(3[2AsmA+ {2~ A^)cosA- 2]/A^- 
-/{-A^ cos A + 4[(3v42 - 6) cos A+ 
{A^-6A) smA + 6])/A'^ 
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where 



(1-^)4 

(l + r?/2)2 



(i = —6rj 
7 = rya/2 



(1 - 



When a second size distribution is included in the fit, 
it has the same C-parametcr and 77 values as the first 
distribution. To set the structure factor to unity, one 
simply sets 77 = 0. This is appropriate for dilute systems. 



D. Fitting routine details 
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FIG. 2: Examples of data from nanoporous silsesquioxane 
films, with different porogen loadings. Open symbols: raw 
data. Lines; fitted curves using SAXSFit. 



The program uses a least-squares fitting routine which 
follows the Levenberg-Marquardt non-linear regression 
method to minimize the reduced x^- The integrals are 
calculated using the Romberg integration method with 
2^" intervals. Since the integral[T]must have finite bounds 
on r, these are chosen based on the range of the distri- 
bution function n(r), such that n{r) < 1 x 10^^^. The 
bounds are calculated numerically as follows: 

For the log-normal and Schulz distributions, the lower 
bound is roexp(— 8cr) and the upper bound roexp(8cr). 

For the Gaussian distribution, the lower bound is the 
maximum of zero or rp — Sti, and the upper bound is 
ro + 8cr. 



C = [J^WJ)-^ 

where J is the Jacobian matrix of partial derivatives and 
is a diagonal matrix where Wn is the weighting on 
the I*'' data point [l^. Finally the reported parameter 
uncertainties are twice the square root of the diagonals 
of C, i.e. 5Pi — 2\JCn. This is two standard deviations, 
which for a Gaussian distribution of errors represents a 
95% confidence interval. 



IV. EXAMPLES 



E. Statistical analysis 

The reduced and (goodness of fit) from the non- 
linear regression are reported at the end of the fitting 
procedure. These are common statistical measures and 
defined as follows: 



Reduced = Wi {y^ - F{xi)f 

n — p ^-^ 

where n is the number of data points, p is the number 
of free parameters, Wi are the weightings, is the input 
data /(g) and F{x) is the calculated I{q). 

{v^~nx)) 



where yi and F(xi) are defined above, and F(x) is the 
average of the F[x) values (a constant). 

The parameter uncertainties are obtained by calculat- 
ing the covariance matrix Cy , from 



Figure 2 shows examples of data and the fitted result 
for nanoporous methyl silsesquioxane films p7| . formed 
by spin-coating a solution of the silsesquioxane along with 
a sacrificial polymer ('porogen'), and then annealing to 
remove the polymer and leave behind a nanoporous net- 
work. As the proportion of porogen is increased, the 
pores are observed to increase in size [17i] . Data are 
shown for films with porogen loadings of 5 to 25 % with 
the background (from methyl silsesquioxane) subtracted. 
A single log-normal size distribution was fitted to each, 
the C-parameter was fixed at 1.1, r] was fixed to the 
porosity (as determined from the porogen loading), and 
no background function was used. The results obtained 
are given in Table |I] and show an increase in the pore size 
with increased porogen loadings, in good agreement with 
electron microscopy and previous results [17j . 

Figure [3] shows data from a nanoporous glass using a 
three arm star shaped polymer as the porogen [l^ , which 
was found to exhibit two pore size distributions. 

The parameters for the fit are as follows: 

First distribution: ro = 50.5 ± 0.5, a = 0.287 ± 0.008. 

Second distribution: rp = 15 ± 1, cr = 0.23 ± 0.05. 

The number fraction of the second distribution was 88 
± 2 %, which equates to a volume of 16 ± 5 %. The 
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Porogen loading 


5% 


10% 


15% 


25% 


Pore size radius (A) 


f8.77 ± 0.f4 


16.94 ± 0.06 


22.42 ± 0.05 


31.67 ± 0.06 


Dispersion 


0.305 ± 0.004 


0.382 ± 0.002 


0.367 ± 0.001 


0.370 ± 0.001 


Reduced 


f.f05 


2.117 


2.913 


2.023 


(degree of fit) 


0.9771 


0.9858 


0.9952 


0.9976 



TABLE I: Fitted parameters for nanoporous silsequioxane film samples shown in Figure [5] The pore radii differ from that 
reported by Huang et al. [l3| due to a slightly different form used for the log-normal distribution. When plotted as a function 
of radius the distributions needed to fit the data are identical. 
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FIG. 3: Example of data from nanoporous glass, showing two 
pore size distributions. Open symbols: raw data. Line: fitted 
curve using SAXSFit. 



(tested on Linux Ubuntu and Mac OS X10.4). The 
SAXSFit programs and user manual are available from 



http: / / www.irl.cri.nz / S AXSfiles.aspx 



VI. SUMMARY 



SAXSFit is a useful program for fitting small angle 
x-ray and neutron scattering data, using a hard sphere 
model with local monodisperse approximation. SAXSFit 
provides an easy-to-use analysis package for novices and 
experts. It is stand-alone software and can be used in a 
variety of software environments. 
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V. SOFTWARE AVAILABILITY AND SYSTEM 
REQUIREMENTS 

SAXSFit is provided as a Windows executable (tested 
on Windows 98, 2000 and XP), or Java .jar executable 
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